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DESCRIPTION 

CODING METHOD AND APPARATUS , AND COMPUTER PROGRAM AND 
COMPUTER -READABLE STORAGE MEDIUM 

5 TECHNICAL FIELD 

The present invention relates to data coding, and 
more particularly, to a coding technique for encoding 
distribution data in multidimensional space such as a 
two-dimensional image and three-dimensional density 
distribution data. 

BACKGROUND ART 

As a general technique for encoding 
multidimensional distribution data such as a two- 
dimensional image, a data amount compression technique 
such as a JPEG compression technique is known. Two- 
dimensional image data is subjected to discrete cosine 
transform then quant ization, and encoding by the 
Huffman coding or the like. 

Further, a method for predictive -encoding image 
data, block- dividing prediction error data by the 
predictive coding, and determining whether or not each 
block has a large amount of prediction error or small 
amount of prediction error is known (for example, 
Japanese Patent Application Laid-open No. Hei 11- 
331852). According to the method disclosed by the 
document, vector quantization is applied to a block 
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with a large amount of prediction error, then the 
difference value between a vector value of a codebook 
retrieved by the vector quantization and the prediction 
error is calculated, then the difference value and the 
5 prediction error in a block with a small amount of 

prediction error are entropy-encoded, and coded data is 
generated by using the entropy code, a block error- 
amount discrimination flag and an index by vector 
quantization. In this manner, from the viewpoint of 
10 space -saving data conservation and high-speed 

communication, a high efficiency data coding technique 

is desired. 

Further, in an arithmetic unit to perform 
discrete convolution on a two-dimensional image 

15 represented by predetermined two-dimensional weight 
distribution data, which is frequently used in image 
processing such as filtering, accelerated computation 
is expected by use of a large-scale parallelized 
product -sum arithmetic elements. However, the product - 

20 sum operation by two-dimensionally parallelized 
arithmetic elements cannot be realized without 
difficulty due to wiring problem and the like . 
Accordingly, a method for decomposing two-dimensional 
weight distribution data into one -dimensional base 

25 group for convolution is desired. 



DISCLOSURE OF INVENTION 
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The present invention has been made in 
consideration of the above situation, and provides a 
high- efficiency data coding technique, particularly for 
encoding distribution data in multidimensional space 
such as a two-dimensional image by decomposing the data 
to several one-dimensional bases. 

According to the present invention, provided is a 
method for encoding n data to be encoded, into data 
representing a base, comprising: 

a setting step of setting the base; 

a reconstruction step of reconstructing n data by 

using the bases; 

an error calculation step of calculating an error 
between the n data to be encoded and the reconstructed 
n data by using a predetermined error evaluation 
function; and 

a correction step of correcting the base based on 

the error. 

Other features and advantages of the present 
invention will be apparent from the following 
description taken in conjunction with the accompanying 
drawings, in which like reference characters designate 
the same name or similar parts throughout the figures 
thereof . 



BRIEF DESCRIPTION OF DRAWINGS 
The accompanying drawings, which are incorporated 
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in and constitute a part of the specif ication, 
illustrate embodiments of the invention and, together 
with the description, serve to explain the principles 
of the invention. 
5 Pig. 1 is a flowchart showing coding processing 

according to a first embodiment; 

Fig. 2 is a flowchart showing coding processing 
according to a second embodiment; 

Fig. 3 is a flowchart showing coding processing 
10 according to a third embodiment; 

Fig. 4 is a flowchart showing coding processing 
according to a fourth embodiment; 

Fig. 5 is a schematic diagram showing an example 
of two-dimensional image including simple rectangles; 
15 Fig. 6 is an explanatory diagram showing a 

conventional structure for transmission of two- 
dimensional weight distribution data; 

Fig. 7 is an explanatory diagram showing a 
structure for transmission of two-dimensional weight 
20 distribution data in a convolution arithmetic unit to 
perform product -sum arithmetic method according to a 
ninth embodiment; 

Fig . 8 is a block diagram showing an apparatus 
construction according to the first embodiment; 
25 Fig. 9 is an explanatory diagram showing a 

structure for transmission of two-dimensional weight 
distribution data in the convolution arithmetic unit to 
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perform the product -sum arithmetic method according to 
the ninth embodiment; 

Fig. 10 shows 8x8 luminance value data 
corresponding to input image in a numerical example 
5 using particular numerical values; 

Fig. 11 shows an initial one -dimensional adaptive 
base group in the numerical example; 

Fig. 12 shows an array of 8x8 numerical values 
corresponding to a reconstructed image from the initial 
10 one-dimensional adaptive base group in the numerical 
example ; 

Fig. 13 shows first correction amounts for a one- 
dimensional adaptive base group with a base number "0" 
in the numerical example; 
15 Fig. 14 shows the one -dimensional adaptive base 

group with the base number "0 M after the first 
correction in the numerical example; 

Fig. 15 shows the one-dimensional adaptive base 
group with the base number "0" after the second 
20 correction in the numerical example; 

Fig. 16 shows the one -dimensional adaptive base 
group with the base number "0 M after the ninth 
correction in the numerical example; 

Fig. 17 is a line graph showing changes of error 
25 E to the corrections on the one -dimensional adaptive 
base group with the base number "0" in the numerical 
example ; 
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Fig. 18 shows the one -dimensional adaptive base 
group with the base number "0" after a convergence of 
the error E in the numerical example; 

Fig. 19 shows the array of 8x8 numerical values 
5 corresponding to the reconstructed image from the 

initial one -dimensional adaptive base group with the 
base number "0", after the convergence of the error E 
in the numerical example; 

Fig. 20 shows an axray of 8x8 numerical values 
10 corresponding to a difference coded image after 
completion of correction of the one- dimensional 
adaptive base group with the base number "0", after the 
convergence of the error E, in the numerical example; 
Fig. 21 shows the first correction amounts for 
15 the one -dimensional adaptive base group with a base 
number "1" in the numerical example; 

Fig. 22 shows the second correction amounts for 
the one-dimensional adaptive base group with the base 
number "1" in the numerical example; 
20 Fig. 23 shows four one -dimensional adaptive base 

groups after completion of correction in the numerical 
example, and arrays of 8x8 numerical values 
corresponding to decoded images reconstructed from the 
base groups; 

25 Fig. 24 shows an example of an input image 

according to the second embodiment; 

Fig. 25 shows an example of decoded images 
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according to the second embodiment; 

Fig. 26 shows an example of decoded images in a 
case where the precision is rounded according to the 
second embodiment; 

Fig. 27 shows an example of a decoded image 
obtained by dividing an input image and performing 
coding and decoding, according to the second 
embodiment ; 

Fig. 28 shows an example of a decoded image 
according to the third embodiment; and 

Fig. 29 shows quantization coefficient vectors 
according to the third embodiment. 

BEST MODE FOR CARRYING OUT THE INVENTION 
15 Preferred embodiments of the present invention 

will now be described in detail in accordance with the 
accompanying drawings . 
<First Embodiment > 

As a first embodiment , a method for encoding a 
20 two-dimensional gray-scale image, defined in a finite 
range, to two one -dimensional adaptive bases will be 
described. Note that in the respective embodiments, 
the adaptive base will also be simply referred to as a 
"base". 

25 Fig. 8 is a block diagram showing an apparatus 

construction according to the first embodiment. 

In Fig., 8, reference numeral 1 denotes a CPU 



5 



10 
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which controls the overall apparatus; 2, a ROM in which 
a processing procedure (program) for the CPU 1 is 
stored; and 3, a RAM used as a work area for the CPU 1, 
in which an image data area 3a for storing image data 
5 to be encoded and an area 3b for storing base data are 
ensured. Note that the sizes of the areas 3a and 3b 
are determined in correspondence with the size of input 
image data. Numeral 4 denotes an image input unit 
which inputs a gray- scale image data to be encoded. As 

10 the image input, an image sensing device or the like 

may be used. In a case where image data is stored in a 
storage medium, a device to access the storage medium 
may be used as the image input unit. Otherwise, in a 
case where an image is received via a network, a 

15 network interface may be used as the image input unit. 
Numeral 5 denotes an external storage unit (hard disk 
or the like) in which coded data is stored. In a case 
where coded data is outputted onto a network, a network 
interface may be used. 

20 In the present embodiment, for the sake of 

simplification of explanation, image data to be encoded 
is handled as a l-pixel-8-bit multivalue (gray-scale) 
monochrome image. (If the image data represents a 
color image, coding to be described below is repeated 

25 in correspondence with the number of components of the 
color image.) Further, the image has a rectangular 
shape where the number of pixels in a horizontal 



WO 2005/074145 



PCT/JP2005/001289 



direction is x max and the number of pixels in a vertical 
direction is y max . In a case where the image input unit 
4 is an image scanner, the number of pixels x ma x and y max 
are uniquely determined respectively based on a reading 
5 resolution and an original size. 

Input image data is stored into the image data 
area 3a. At this time, the pixel value of a coordinate 
(x,y) (Note that Ossx<x max and 0^y<y max hold) is expressed 
as f(x,y). For example, in a case where data in k- 

10 dimensional space is handled, the data is expressed as 
f(x x , x 2 , x k ). Further, as the size of the input 

image data is determined, two one-dimensional adaptive 
bases X(x) and Y(y) are allocated in the base area 3b. 
Each one -dimensional adaptive base can be regarded as a 

15 continuous or discrete one -dimensional sequence of 

numerical values. Further, if the continuous sequence, 
of numerical values are amplitude values, the sequence 
of numerical values can be regarded as a waveform. The 
ranges of n x n and "y" are the same as those of the 

20 input image. Note that the pixel value f (x,y) is an 8- 
bit (1 byte) integer value, and the one -dimensional 
adaptive bases X(x) and Y(y) are ensured in multi-byte 
format in the base area 3b of the RAM 3 for storage of 
decimal place values. 

25 Hereinbelow, a coding processing procedure of the 

CPU 1 according to the first embodiment will be 
described with reference to the flowchart of Fig. 1 
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( the corresponding program is stored in the ROM 2 ) • 

First, in image data input (step S10), an image 
f (x,y) as a two-dimensional gray-scale image to be 
encoded is inputted and stored in the image area 3a. 
5 Next, base setting is performed (step Sll). At this 
step, initial values of the two one-dimensional 
adaptive bases X(x) and Y(y) are set. Note that as the 
one -dimensional adaptive bases X(x) and Y(y) , values 
which are not "0" in the limited ranges 0^x<x max and 

10 0<sy<y max are set. In this initial stage, the values are 
set by using random functions. Note that in this 
embodiment, the respective definition ranges are 
limited as above, however, the definition ranges are 
not limited to such limited ranges. That is, even if 

15 f(x,y), X(x) and Y(y) are defined in indefinite 

definition ranges, the present method is applicable 
(however, there are limitations since the memory 
capacity is finite) unless an error E and correction 
amounts AX(x) and AY(y) to be described later diverge. 

20 In base reconstruction processing (reconstruction 

processing from a group of adaptive base) (step S12) , a 
reconstructed image f(x,y) is generated (reconstructed 
image from a group of adaptive bases) as follows. 

25 f(x,y)=x(*>Y(y) 

Next, in error calculation processing (step S13) , 
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an error E between the input image f (x,y) and the 
reconstructed image f(x,y) is calculated by using the 
following error evaluation function. 



Note that /f (x)dx means an integrated value of 
f(x) in the entire x definition range. As integration 
is similarly performed regarding the y component, the 

10 error E is an accumulated error in the entire area of 
the input image. Further, as the error E is used in 
error convergence determination (step S16 as branch 
determination) to be described later, it is 
sequentially stored as a record in the RAM 3. 

15 Next, based on the error E, correction amounts 

for the respective one-dimensional adaptive bases are 
calculated (step S14). The correction amounts AX(x) 
and AY(y) are calculated as follows. 



5 




20 




Note that correction parameters v\ x and r\ Y are 



positive constants (setting of the parameters r| x and t) y 



will be described later). 
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The functions to obtain the correction amounts 
are partial differentiation of the expression to obtain 
the error E . That is , the error E always has a value 
greater than "0", and in the ultimate sense, the object 
5 of the present embodiment is to obtain one -dimensional 
adaptive bases with a minimum error E. In a case where 
the slope in the current one-dimensional adaptive base 
is minus, a minimum value one-dimensional base is 
positioned in the positive direction from the current 

10 point, while in a case where the slope in the one- 
dimensional adaptive base is plus, a minimum value one- 
dimensional base is positioned in the negative 
direction form the current point. Accordingly, 
multiplying the slope by minus gives a correction 

15 amount corresponding to the slope. 

Next, at step S15, correction processing is 
performed on the respective bases. That is, the 
correction amounts AX(x) and AY(y) obtained at step S14 
are added to the uncorrected one -dimensional adaptive 

20 bases X(x) and Y(y) , thereby corrected one -dimensional 
adaptive bases X'(x) and Y'(y) are obtained. The 
processing is expressed as follows. 
X 1 (x)=X(x)+AX(x) 
Y f (y)=Y(y)+AY(y) 

25 Upon correction of one -dimensional adaptive bases, 

the above-described condition "the respective one- 
dimensional adaptive bases are always not "0" as long 
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as the Input image f (x,y) is always not "0" in the 
entire defined x and y ranges" must be satisfied. Note 
that if such condition is satisfied, minute 
fluctuations can be applied to at least a part of the 
5 one-dimensional adaptive bases in a random manner. 

The above -described correction on one -dimensional 
adaptive bases at step S12 (base reconstruction 
processing) to step S15 (base correction processing) is 
repeatedly performed until completion determination is 
10 made at step S16 at which error lower limit is 

determined and at step 17 at which error convergence is 
determined. 

In the error lower limit determination at step 
S16, if the error E of immediately previous correction, 

15 calculated in the error calculation (step S13), is 
equal to or less than a predetermined value, it is 
determined whether or not the one-dimensional adaptive 
bases produce an error E within the allowable range. 
If the above completion determination is made, the one- 

20 dimensional adaptive bases at that point are recorded 
as coded data, and the process ends. 

Further, in the error convergence determination 
at step SI 7, it is determined whether or not the error 
has converged based on the recorded error E calculated 

25 in the error calculation (step S13). If it is 
determined that the error has converged, the 
determination that the one -dimensional adaptive bases 
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at that time produce a minimum error value is outputted, 
and the process ends. The one-dimensional adaptive 
bases upon determination of convergence are recorded as 
coded data of the input image f(x,y). As the 
5 determination of error convergence, if a status where 
the error E has almost no change (the error difference 
is equal to or less than a predetermined threshold 
value) in each correction has continued, it is 
determined that the error has converged. 

10 By the above -de scribed processing, the two- 

dimensional multivalue (gray-scale) image f (x,y) can be 
encoded to data expressed by two one -dimensional 
adaptive bases X(x) and Y(y) . In other words, as the 
image data f (x,y) having pixels x max xy max is reduced to 

15 the coded data expressed by Xmax+Ymax. it can be assumed 
that the initial data is compressed (note that since 1 
pixel of the input image corresponds to 8 bits, while 1 
element of one -dimensional adaptive base corresponds to 
several bytes, the compression ratio cannot be 

20 specified simply by the number of data) . Further, a 

decoded image f decode (x,y) corresponding to the coded two 
one-dimensional adaptive bases X(x) and Y(y) is 
obtained based on the following expression. 
fdecode(x,y)=X(x)xY(y) 

25 In the above coding processing, the initial one- 

dimensional adaptive bases X(x) and Y(y) set in the 
base setting at step Sll are, preferably, values about 
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within the defined x and y ranges of the input image. 
Preferably, if values which are not 11 0" but 
approximately "0" are set, initial value dependency is 
reduced and appropriate one -dimensional adaptive bases 
5 can be obtained with high probability. 

Further, preferably, the correction parameters Y) x 
and Tiy are determined in correspondence with the range 
of the input image f(x,y) and determined x and y ranges. 
More particularly, it is preferable that the parameter 

10 Y| x is in inverse proportion to the second power of the 
y range and the range of the input image f(x,y), and 
the parameter r\ y is in inverse proportion to the second 
power of the x range and the range of the input image 
f(x,y). In this embodiment, the correction parameters 

15 Yj x and r\ Y are changed in accordance with the respective 
one -dimensional adaptive bases, however, the same value 
may be used if the defined x and y ranges are 
approximately the same. If the correction parameter is 
set to a great value, the completion determination in 

20 the error lower limit determination (step S16) or in 
the error convergence determination (step S17) can be 
made with high probability by a reduced number of 
correction of one -dimensional adaptive bases. However, 
there is a high probability that the error E greatly 

25 fluctuates around a minimum value and appropriate one- 
dimensional adaptive bases cannot be obtained. Further, 
if the correction parameter is set to a small value. 
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the error E around a minimum value can be obtained and 
appropriate one -dimensional adaptive bases can be 
obtained with high probability. However, there is a 
high probability that a large number of correction of 
5 one-dimensional adaptive bases is required until the 
completion determination is made in the error lower 
limit determination (step S16) or in the error 
convergence determination (step S17). In view of the 
above points, it is preferable that the correction 

10 parameter is set to a small value not to cause great 
fluctuation of the error E around a minimum value and 
not to require an extremely large number of correction 
of one- dimensional adaptive bases. Further, a method 
to gradually reduce the correction parameter in 

15 accordance with the number of corrections may be 
applied. 

Further, it is preferable that the previously- 
determined allowable range of the error E, used in the 
error lower limit determination (step S16), is 

20 determined in correspondence with the range of the 

input image f(x,y) and the defined x and y ranges. If 
the error E is set to a great value, the completion 
determination in the error lower limit determination 
(step S16) can be made with high probability by a 

25 reduced number of correction of one -dimensional 

adaptive bases. However, there is a high probability 
that the one -dimensional adaptive bases with high error 
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E value are obtained- On the other hand. If the error 
E ±s set to a small value, appropriate one -dimensional 
adaptive bases with lower error E value are obtained 
with high probability. However, there is a high 
5 probability that a greater number of correction of one- 
dimensional adaptive bases is required by the 
completion determination in the error lower limit 
determination (step S16) or completion determination is 
not made in the error convergence determination (step 
10 S17). 

Further, in a case where the correction amounts 
AX(x) and AY(y) are calculated in the above- described 
correction amount calculation (step S14) as follows, by 
using previous correction amounts AX(x),AY(y) , there is 
15 a high probability that more appropriate one- 

dimensional adaptive bases with less error value cab be 
obtained, or desired one -dimensional adaptive bases can 
be obtained with a reduced number of corrections. 

A*to~»7* f]f(x,y)-f(pc,y)\'Y(y)dy +cc x -&X{x) 
AF(y)-»7 y f\f(x,y)-f(x,y)\x(x)dx + a y -Ar(y) 

Note that momentum constants a* and oty are 
integer constants less than "1". In the present 
embodiment, the momentum constants are changed in 
accordance with the respective one -dimensional adaptive 
bases, however, the same value may be used. If the 



20 



25 
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momentum constant Is set to a great value, the above 
advantages are obtained with high probability, however, 
adverse effects also appear and the correction of one- 
dimensional adaptive bases may become unstable with 
5 high probability. Accordingly, the momentum constant 
set to an extremely large value is not preferable. On 
the other hand, if the momentum constant is set to a 
small value, the probability of appearance of adverse 
effects is lowered, however, the probability to obtain 

10 the above advantages is lowered. Accordingly, it is 

preferable that the momentum constant is set to a value 
as large as possible within a range where the 
correction of one-dimensional adaptive bases is stable. 
The first embodiment has been described as an 

15 example of coding of two-dimensional gray- scale image 

f (x,y) into two one-dimensional adaptive bases X(x) and 
Y(y) , however, the image to be encoded is not limited 
to a two-dimensional image. For example, the above 
coding is applicable to density distribution in three- 

20 dimensional space, or applicable to multidimensional 
space higher than the three-dimensional space. In a 
case where an object in three-dimensional space is 
encoded, the object can be encoded to three one- 
dimensional adaptive bases X(x) , Y(y) and Z(z). 

25 

< Second Embodiment > 

In the above-described first embodiment, as long 
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as a two -dimensional image to be encoded is 
comparatively simple, error-reduced coding can be 
realized by appropriately selecting a coordinate system. 
That is, as long as the two-dimensional image comprises 
5 simple rectangular shapes as shown in Pig. 5, coding 

with comparatively small error can be realized by using 
an orthogonal coordinate system orthogonal to the 
rectangles. In this manner, by appropriately selecting 
a polar coordinate system or the like in correspondence 

10 with data, or by selecting a columnar coordinate system 
or the like for three-dimensional data, coding with 
reduced error E can be realized. 

However, in a case where a two-dimensional image 
which is a natural image is encoded by the method 

15 according to the above-described first embodiment, the 
error between data (image data) decoded by using coded 
data and original data becomes large. Accordingly, in 
the present embodiment, a method for encoding a 
discrete two-dimensional gray-scale image f (x,y) 

20 defined in a finite range into plural group of discrete 
one -dimensional adaptive base groups Xi(x), Yi(y)(x and 
y are integers) will be described as a modification to 
the first embodiment. 

The outline of the second embodiment will be 

25 described. One one-dimensional adaptive base group is 
subjected to similar processing to that of the first 
embodiment, and subjected to the lower limit 
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determination of error E and convergence determination. 
Then, if it is not determined that the error E is 
sufficiently small , a new one -dimensional adaptive base 
group is added (allocated), and the same processing is 
5 repeated so as to reduce the error between the coded 
data and decoded data. 

In the present embodiment, the apparatus 
construction is the same as that of the first 
embodiment. Hereinbelow, the processing according to 
10 the second embodiment will be described with reference 
to the flowchart of Fig. 2. 

First, in image input (step S20) , an input image 
f(x,y) as a discrete two-dimensional gray-scale image 
to be encoded is inputted. As in the case of the first 
15 embodiment, x and y have predetermined definition areas, 
0£x<x raa x and 0sy<y max - Further, in the second embodiment, 
the input image is discrete, and x and y are integers. 

As a particular example, processing performed on 
8x8 luminance value data as shown in Fig. 10 will be 
20 described. Hereinbelow, the example using particular 
numerical values will be referred to as a numerical 
example. The numerical values in Fig. 10 are non- 
negative integers in 8 -bit precision, however, any data 
may be employed as long as difference between the 
25 values can be defined. Following the image input (step 

520) , the base setting processing is performed (step 

521) . One of discrete one-dimensional adaptive base 
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group Xi(x), Yi(y)(x and y are integers) is selected 
and set as an initial one -dimensional adaptive base 
group. The suffix "1" is a base number. In the 
present embodiment, although the base number is not 
5 particularly limited, the base number is sequentially 
selected from "0" to a predetermined number. As the 
image subjected to processing in the numerical example 
is an 8x8 luminance value data, each of the one- 
dimensional adaptive bases is a sequence of numerical 

10 values having 8 elements. As described in the first 
embodiment, it is preferable that approximately "0" 
random numbers are set as initial numerical values. 
However, in the initial one -dimensional adaptive base 
group having suffix "1". a one-dimensional adaptive 

15 base group as shown in Fig. 11 is used merely to 

satisfy the condition that all the elements are not "0" , 
for the sake of simplification. In Fig. 11, the value 
of bracketed "t" indicates the number of corrections on 
one -dimensional adaptive base. That is, if t=n holds, 

20 the one -dimensional adaptive base has been corrected n 
times. 

Next, the base reconstruction processing is 
performed (step S22). A reconstructed image f x (x,y) is 
calculated as follows. 

25 
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Fig. 12 shows an array of 8x8 numerical values 
corresponding to a reconstructed image f Q (x,y) 
immediately after the base setting processing (step 
S21) using the initial one -dimensional adaptive base 
group. Next, the error E is calculated as follows by 
error calculation (step S23) . 



(1) 



10 In the above expression, fi(x) is a difference 

coded image expressed as follows. 

f,(.x,y)-f(x,y)-y~f„(x,y) 

15 Note that ]>)/(*) means summation of f (x) when x 

X 

is an integer, in the entire x definition range. 

m 

Further, means liner summation of fk(x) 

when n sdcsan holds. If m<n holds, the result of 
summation is 0. In the numerical example, as the 
20 difference coded image f 0 (x) becomes an input image, 
the error E is "482913" . 

The above difference coded image fi(x,y) means 
the result of subtraction of the respective 
reconstructed images, reconstructed from one one- 
25 dimensional adaptive base group, having the base number 
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"0" to "1-1", from the input image f(x,y). That is, in 
the second embodiment, as the respective one- 
dimensional adaptive base groups are corrected while a 
one-dimensional adaptive base group is sequentially 
5 selected from the base number n 0'\ the difference coded 
image fi(x,y) corresponds to the result of subtraction 
of the respective reconstructed images, reconstructed 
from the correction -completed one -dimensional adaptive 
base groups, from the input image f(x,y). 
10 Next, in the correction amount calculation (step 

S24), correction amounts AXi(x) and AY x (y) for the one- 
dimensional adaptive base group X x (x) . Y x (y) , selected 
in the base setting (step S21), are calculated as 
follows . 

15 

a*, to- -n, -^j - n,-2fo(*,y)- f.Mfrb)] 
at to- -n, - v, 2fob,y)-f,kyi *,to] 

Otherwise, the correction amounts may be 
calculated by using the following expressions where the 

A A 

20 previous correction amounts AX(x),AY(y) are multiplied 
by the momentum constants a*, cty described in the first 
embodiment are added. 
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. . -(2) 

In the above expressions, Tlx and riy are 
5 correction parameters similar to those described in the 
first embodiment. 

In the numerical example , the correction amounts 
are calculated by using the expression (2). In the 
case of correction amounts for the one- dimensional 
10 adaptive base group shown in Fig. 11, i.e., immediately 
after the base setting processing (step S21), the 
values of the momentum constants are "O" since all the 

A A 

elements of the previous correction amounts AX(x),AY(y) 
are "0". Fig. 13 shows the correction amounts in this 

15 case. In the figure, the value of bracketed M t n 

indicates the number of corrections. That is, if t=n 
holds, the correction amounts corresponds to the n-th 
one -dimensional adaptive base correction. As the 
correction parameters in this case, Ti x =T) y =0 . 00078 holds, 

20 and as the momentum constants, ctx=cxy=0 . 4 holds. 

In the base correction (step S25), assuming that 
corrected one -dimensional adaptive base group is X f x (x), 
Y , i(y), and uncorrected one-dimensional adaptive base 
group is X x (x) , Y x (y) , the respective one -dimensional 
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adaptive base groups are corrected by using the 
correction amounts AXi(x) and AYi(y) calculated as 
above as follows ♦ 

A?,(r) 

5 y,ty)«^(y)+Ar,(y) 

The result of this correction on the one- 
dimensional adaptive base group shown in Fig. 11 is as 
shown in Fig. 14. In the first correction, since all 

10 the elements of the initial one -dimensional adaptive 
base group are "0" except one element, only a part of 
the elements has been corrected. Further, as the 
numerical value other than n 0" is comparatively small, 
"0.1", the first correction amounts are very small. 

15 Preferably, random numbers are set as the numerical 

values, and all the elements are corrected in the first 
correction. However, since it is desirable that 
initial values are approximately "0", the correction 
amounts are small. Then the error E is calculated 

20 based on the expression (1) by using the corrected one- 
dimensional adaptive base group, and the value 
"482911.48" is obtained. The error E is slightly 
reduced in the first correction. As in the case of the 
first embodiment, such correction is repeatedly 

25 performed until the completion determination is made in 
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the error lower limit determination (step S26) or the 
convergence determination is made in the error 
convergence determination (step S27). thereby the one- 
dimensional adaptive base group having the base number 
selected in the base setting (step S21) is corrected. 
For example, after the correction with the first 
correction amounts as shown in Fig. 13, the second 
correction amounts are as shown in Fig. 15, and the 
error E is "482906.05". In the second correction, 
although the correction amounts are very small, all the 
elements are corrected. In this numerical example, the 
seventh to tenth correction amounts are comparatively 
large. For example, the ninth correction amounts are 
as shown in Fig. 16. The error E immediately before 
the ninth correction is "286696.02", and the error 
immediately after the ninth correction is greatly 
reduced to "16349.27". Fig. 17 shows the relation 
between the corrections and the changes of the error E. 
As shown in Fig. 17, the error E is reduced by each 
correction in an accelerative manner, from the first to 
ninth corrections. Then the error is temporarily 
increased, then, the error, with fluctuations, is 
gradually converged. The sudden increase after the 
ninth correction is caused by the influence of momentum 
constants. If the momentum constants are set to "0", 
such increase does not easily occur, however, above 
processing is effective for reduction of the number of 
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corrections by the convergence and for avoidance of a 
trap to a local minimum value. 

The error lower limit determination at step S26 
and the error convergence determination at step S27 are 
5 the same as those described in the first embodiment. 
In the error lower limit determination at step S26 , if 
completion determination has been made, the number of 
base groups selected and corrected by that time and 
those one-dimensional adaptive base groups are recorded 

10 as coded data, and the process ends. In the above 
numerical example, the error is still large at this 
stage and sufficient approximation is not made. For 
example, if the error lower limit is "100", the 
completion determination is not made. If the 

15 completion determination has not been made, the process 
proceeds to the next error convergence determination at 
step S27, at which it is determined whether or not the 
error E has converged. In the numerical example, when 
correction has been performed about 30 times, the error 

20 E has approximately converged, to about "16199". Fig. 
18 shows the one -dimensional adaptive base group when 
the error E has converged. Further, Fig. 19 shows the 
array of 8x8 numerical values corresponding to the 
reconstructed image f Q (x,y) from the one -dimensional 

25 adaptive base group. In the figure, the number of 
corrections of one-dimensional adaptive base "t-oo" 
means that correction has been repeated such that the 
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error E has converged. 

In the first embodiment, if it is determined that 
the error E has converged in the error convergence 
determination at step S27, completion determination is 
5 made and the process ends. However # in the second 
embodiment, the process proceeds to difference 
calculation processing at step S28. If it is not 
determined that the error E has converged, the process 
returns to the base reconstruction at step S22, in 
10 which correction is repeatedly performed on the one- 
dimensional adaptive base group with the selected base 
number . 

If it is determined that the error E has 
converged, then in the difference calculation at step 
15 S28, a reconstructed image f^y) reconstructed from 

the one -dimensional adaptive base group with the base 
number "l n is calculated by 



20 

and a difference coded image fi+i(x,y) is generated by 





25 In the above -de scribed numerical example, as 1=0 

holds, the difference coded image fi(x,y) is the result 
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of simple subtraction of the reconstruction image 
f 0 (x,y) , reconstructed from the one-dimensional adaptive 
base group with the base number "0" , from the input 
image f (x # y) . Fig. 20 shows an array of 8x8 numerical 
5 values corresponding to the difference coded image 
fi(x,y) - 

Next, in the base number determination at step 
S28, if it is determined that the number of correction 
completed base groups has become a predetermined number , 

10 completion determination is made. Then the number of 
selected and corrected base groups by that time and 
those one -dimensional adaptive base groups are recorded 
as coded data, and the process ends. If it is 
determined that the number of correction completed base 

15 groups has not become the predetermined number, the 
process returns to the base setting at step S21, in 
which a new one-dimensional adaptive base group is 
selected, and the processing from the base 
reconstruction at step S22 to the error convergence 

20 determination at step S27 is repeatedly performed, thus 
correction on the newly- selected one -dimensional 
adaptive base group is performed. 

In the above -described numerical example, since 
correction has been completed regarding only one group, 

25 if an allowable number of base groups is set to, e.g., 
"4", as the number of correction completed groups has 
not become the predetermined number, the process 
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returns to the base setting (step S21). Then the next 
one -dimensional adaptive base group having the base 
number "1" is selected, and correction is performed on 
the one -dimensional adaptive base group in accordance 
5 with the expression (2) by using the difference coded 
image f x (x,y) obtained in the difference calculation 
(step S28) . 

As described above, as the initial one- 
dimensional adaptive base group is as shown in Fig. 11, 

10 the one -dimensional adaptive base group with the base 
number "1" selected here is the same as that shown in 
Fig. 11. At this time, the first correction amounts 
are as shown in Fig. 21, and only a part of the 
elements is corrected as in the case of the first 

15 correction on the one- dimensional adaptive base group 
with the base number "0". The immediately previous 
error E of the first correction is "16199 . 19" , and the 
error E becomes "16199.14" by the correction, thus the 
error is slightly reduced. The second correction 

20 amounts are as shown in Fig. 22. The error E after the 
correction is "16199.07". Such correction is repeated 
by the error convergence, thereby the error E is 
reduced to about "4749". 

When the error has converged, the next difference 

25 coded image f 2 (x,y) is obtained based on the expression 
(3). The difference coded image f 2 (x,y) is the result 
of subtraction of the reconstruction image f 0 (x,y) 
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reconstructed from the one-dimensional adaptive base 
group with the base number M 0" and the reconstruction 
image f x {x 9 y) reconstructed from the one -dimensional 
adaptive base group with the base number n l" f from the 
5 input image f(x,y). 

In this manner, the difference coded image is the 
data component obtained by subtracting the 
reconstruction images reconstructed from all the 
correction- completed one -dimensional adaptive base 
10 groups from the input image. With respect to the 
difference coded image f 2 (x,y) * the initial one- 
dimensional adaptive base group with the base number 
"2" is selected, and correction is performed on the 
one-dimensional adaptive base group based on the 
15 expression (2). In this manner, the number of one- 
dimensional adaptive base groups is increased. 

In the numerical example, when the correction of 
the first one -dimensional adaptive base group has been 
completed, the error E is about " 16199" , and when the 
20 correction on the second base group has been completed, 
the error E is about n 4749 M . Then, when the correction 
on the third base group has been completed, the error E 
is about w 3312" , and when the correction on the fourth 
base group has been completed, the error E is about 
25 M 671" . In this numerical example, as the allowable 

number of the base groups is "4", the process ends at 
this stage by the determination of the number of base 
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groups (step S28) regardless of the error E value. In 
a case were the error lower limit in the error lower 
limit determination (step S26) is set to "1000" . as the 
error becomes less than the set lower limit in the 
5 middle of the fourth correction, the process ends at 
that time. 

However, the above arrangement, where the 
correction is terminated in the middle without 
increasing the number of base groups even though the 

10 error can be further reduced, is significant in the 

point of high-speed processing, but not advantageous in 
the point of approximation performance. Accordingly, 
the completion determination in the error lower limit 
determination (step S26) may be arranged such that, in 

15 correction on one -dimensional adaptive base group where 
the error becomes less than the lower limit, completion 
determination is made when the error has converged. 

The predetermined number of base groups used in 
the number of base groups determination processing 

20 (step S28) may be set to an infinite number, however, 
if the error allowable range used in the error lower 
limit determination (step S26) is small, there is a 
high probability that a very large number of base 
groups are required. Accordingly,- the number of base 

25 groups is previously limited, thereby the above coding 
can be performed with data amount with upper limit. 

The coded data, encoded by the above method, is 
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decoded to a decoded Image f decode (x,y) as follows. 

In the above expression, "k M means the number of 
5 selected and corrected base groups. 

In this manner, in the coding method according to 
the second embodiment, decoded data can be obtained 
from coded data encoded by simple processing. 

Fig. 23 shows four one -dimensional adaptive base 
groups obtained in the numerical example, and arrays of 
8x8 numerical values corresponding to decoded images 
reconstructed from the base groups. 

In the numerical example, the image to be encoded 
is an 8x8 luminance value image for the sake of 
simplification, however, such small data size has 
almost no meaning in the point of data compression. It 
is preferable that coding is performed on larger sized 
data as an example of image data coding for this 
purpose . 

For example. Fig. 25 shows reconstructed images 
obtained by performing the above coding on a 100x100 
pixel image as shown in Fig. 24 and reconstructing from 
the coded image. 

In Fig. 25, the reconstructed images are obtained 
from processings using one base group, two base groups, 
five base groups, ten base groups, and twenty base 
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groups. As it is understood from this example, 
sufficient approximation to the original input image 
can be attained in the case of about twenty base groups. 
If one pixel value of the 100x100 pixel image as data 
5 to be encoded is a 4-byte floating point, the data 

amount is 4 bytesxl00xl00=40000 bytes. If each element 
value of one -dimensional adaptive base is also a 4-byte 
floating point, the data amount of twenty base groups 
is 4 bytesx(100[the number of elements of X-directional 
one-dimensional adaptive bases 3+100 [the number of 
elements of Y-directional one -dimensional adaptive 
bases] )x20( the number of base groups) =16000 bytes. 
Thus, the data compression has been achieved. 

Actually, in the case of general image as shown 
in Fig. 24, as one pixel value is represented by 1 byte 
the data amount is 1 bytesxl00xl00=10000 bytes. In 
this case. If each element value of one -dimensional 
adaptive base is also a 4-byte floating point, the data 
amount is increased. However, in the case of general 
image, as the precision of 4 bytes is not necessary for 
the value of each element of one -dimensional adaptive 
base, precision of 8 bits or 6 bits is sufficient. 
Accordingly, it may be arranged such that only upper 
and lower limits of the respective one -dimensional 
adaptive bases are held with the 4-byte floating point 
precision, and the values of elements of the one- 
dimensional adaptive base are rounded to a desired 
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number of bits by using the upper and lower limits. 

Fig. 26 shows decoded images reconstructed from 
coded images in a case where the number of base groups 
is twenty and the bit numbers of adaptive bases are 
5 respectively rounded to 8-bits, 6-bits and 4-bits. The 
data amount of the rounded 8 -bit image is {1 byte [8 
bits ]x( 100 [the number of X-directional one -dimensional 
adaptive base elements] +100 [the number of Y-directional 
one-dimensional adaptive base elements ]+4bytesx 4 [upper 

10 limit necessary for 1 base group/number of held lower 
limit data] }x20[the number of base groups]=4320 bytes. 
Thus, the data amount is compressed. In the case of 6- 
bit image, the data amount is 3320 bytes, and in the 
case of 4-bit image, the data amount is 2320 bytes. 

15 The rounding of the values of one -dimensional 

adaptive base elements is not limited to the above - 
described method, but any other method such as rounding 
using a predetermined absolute value may be employed. 
Further, the rounding of the values of one -dimensional 

20 adaptive base elements may be performed after the 
completion of coding using plural one -dimensional 
adaptive base groups, however, there is a probability 
that the error is accumulated. Accordingly, it is 
preferable that the rounding is performed upon 

25 completion of correction on one one -dimensional 

adaptive base group then a difference coded image to be 
used for correction on the next one -dimensional 
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one -dimensional base group. 

Further, in a case where the above -described 
coding method is applied to a larger sized image of, 

5 e.g., 1280x960 pixels, as the number of one-dimensional 
adaptive base elements is large, the error convergence 
cannot be easily attained. In such case, it may be 
arranged such that the image is divided into blocks 
having a predetermined size of, e.g., 50x50 pixels, and 

10 the above- described coding is performed on the 

respective blocks. Fig. 27 shows a decoded image 
obtained by performing the 50x50 pixel size block 
division and coding on the image in Fig. 24 and 
decoding on the coded image. In this case, in an 

15 upper-left 50x50 pixel block, the number of base groups 
is 8, and in upper-right, lower-left and lower-right 
blocks, the numbers of base groups are 10, 7 and 12. 
The one -dimensional adaptive base elements are 
respectively rounded to 8 bits. The size of coded data 

20 is 4292 bytes including upper-limit and lower-limit 
retention data. In comparison with the case of 
encoding a 100x100 sized image with 8-bit rounding 
without division, the data amount is somewhat smaller. 
In this manner, as an input image is divided into 

25 appropriate sized blocks and coding is performed on the 
respective blocks, coding with a small number of base 
groups can be realized in a simple structured block. 
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and the data amount can be reduced. Note that if the 
block size is too small, the advantage of data 
compression by each block is reduced as in the case of 
the above -described 8x8 image, and further, the upper 
5 and lower limit data are increased. Accordingly, it is 
preferable to set an appropriate block size. 

The block size may be adaptively determined for 
an image to be encoded, but the block size may be 
previously set to a predetermined size. In the case 

10 where the block size is a predetermined size, 

remainder (s) may occur upon division. For example, in 
a case where an input image having a size of 1280x960 
pixels is divided by 50x50 pixels from an upper left 
corner, horizontal 30 pixels remain in a right end 

15 portion and 10 pixels remain in a lower end portion. 

In this case, coding may be performed on the right end 
portion on the assumption that it is divided by 30x50 
pixels, on the lower end portion on the assumption that 
it is divided by 50x10 pixels, and on a lower right 

20 portion on the assumption that it is divided by 30x10 
pixels. Further, in the case of coding on all the 
image divided by 50x50 pixels, coding may be performed 
on the assumption that 40 pixels in the lower end 
portion, 20 pixels from the left and 40 pixels from the 

25 bottom in the lower right portion are "0". 

As described above, in the second embodiment, a 
discrete two-dimensional gray- scale image f (x,y) 
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defined in a finite range is encoded to plural discrete 
one-dimensional adaptive base groups Xi(x), Yi(y). As a 
result, in comparison with the first embodiment, even 
in a case where an image to be encoded is a natural 
5 image, coding with high precision can be performed. 

< Third Embodiment > 

In a third embodiment as a modification to the 
above -described second embodiment, a discrete two- 
10 dimensional gray-scale image f (x,y) defined in a finite 
range is divided into blocks and converted to plural 
discrete one- dimensional adaptive base groups X^i(x) , 
^i(y) by each block, and coding is performed on each 
block. In this embodiment, as the apparatus 
15 construction is the same as that of the first 

embodiment, the explanation thereof will be omitted. 

Fig. 3 is a flowchart showing the processing 
according to the third embodiment. Hereinbelow, a 
processing procedure according to the third embodiment 
20 will be described with reference to Fig. 3. 

First, in the image input (step S30), an input 
image f (x,y) which is a discrete two-dimensional gray- 
scale image to be encoded is inputted. As in the case 
of the first and second embodiments, the definition 
25 areas of x and y are integers having ranges 0^x<x max and 

0iSy<ymax. 

Next, in image division processing (step S31), 
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the input image is divided into plural windows W^acy) 
having a size of L x , L y (L x and L y are non-negative 
integers), expressed as follows. 

W**fcy)-f 

Note that i and t> are integers equal to or 
greater than n 0" . The definition ranges of each window 
is Ossx<L x , 02sy<L y . Further, in a window where x+^L x srx max 
or y+tLy^Ymax holds, it may be arranged such that the 
range of x+£L x ^x max or y+£L y sy ma x is rounded down as a 
window having a size smaller than L x , L y . 

Next, window selection processing is performed 
(step S32). One window W*' 5 (x,y) is selected from the 
windows obtained by dividing the input image by the 
image division processing at step S31. The order of 
selection is not particularly limited as long as one 
un- selected window is selected. 

Next, base grouping processing (step S33) is 
performed on the window W^(x,y) selected in the window 
selection processing at step S32. The selected window 
W^' e (x,y) is converted into plural discrete one- 
dimensional adaptive base groups X i ' t 1 (x) , Y % \(y) • Th © 
base grouping processing may be performed by using a 
similar method to that of coding to plural one- 
dimensional adaptive base groups as described in the 
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second embodiment. The definition ranges of one- 
dimensional adaptive base groups X l ' e i(x) , Y*' C i(y) are 
basically 0ssx<L x and 0<y<L y . Further, in the above- 
described window where x+^Lx^Xmax or y+^Ly&Vmax holds, it 

5 may be arranged such that, assuming that the window 
u u u u 

size is L x ,L y , Q*x<L x and 0ay<L y hold. The image 

division processing at step S31 is, more particularly, 

similar to the division by a predetermined block size 

described in the second embodiment. 

10 Since the base grouping processing at step S33 is 

similar to that of the second embodiment, the detailed 
explanation thereof will be omitted. When the 
conversion of the one window vfi x (x,y) , selected in the 
window selection at step S32 to one -dimensional 

15 adaptive base groups by the base grouping processing at 
step S33, has been completed, the process proceeds to 
selection completion determination processing at step 
S34. 

In the selection completion determination 
20 processing at step S34, it is determined whether or not 
all the windows have been selected. If there is an 
unselected window, the process returns to the window 
selection at step S32, in which the unselected window 
is selected, then the conversion of the window to 
25 plural one -dimensional adaptive base groups is 

performed by the above- described base grouping at step 
S33. If it is determined all the windows have been 
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selected, the process proceeds to DCT transform 
processing at step S35. 

In the DCT transform processing at step S35, each 
of the plural discrete one -dimensional adaptive base 
5 groups X^\(x) , Y^i(Y) of each window W^(x,y) , 

generated in the base grouping processing at step S33, 
is discrete-cosine transformed to DCT vectors U^i(x), 
V^i(y). The discrete cosine transform on the 
respective one -dimensional adaptive bases is performed 
10 by using the following transform expression. 

When x=0 holds, ^=1/^2- Further, when x*0 holds, 
15 c x =l holds. Similarly, when y=0 holds, c y =1/^2 holds, 
and when y?*0 holds, c y =l holds. 

In the above transform, all the discrete one- 
dimensional adaptive base groups X ?, S.(x) , ^\(y) 
generated in the base grouping processing at step S33 
20 are transformed to DCT vectors U l,e a(x) , V^i(y). 

Next, in quantization processing at step S36, 
quantization is performed on all the DCT vectors 
U^'^i(x) , V*\(.y) obtained in the DCT transform 
processing at step S35, thereby quantized bases U q l ' t i (x) , 
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V q l '^i(y) are generated. The quantization is performed 
as follows by using one -dimensional quantization 
coefficient vectors Q«(x) and Q y (y) having the same 
number of DCT vector elements. 



U s *,(x)\ 

V"i{y)\ 



|#J means a maximum integer not greater than a. 
For example, [15J=1 holds. In the present embodiment, 
rounding down is performed as above, however, the 
integer may be obtained by rounding off. 

As in the case of JPEG quantization coefficients, 
it is preferable that the quantization coefficient 
vectors Q*(x) and Q y (y) are increased as the x 

15 component is increased. Further, the values of the 

quantization coefficient vectors Q x (x) and Q y (y) may be 
previously determined in correspondence with the ranges 
of DCT vectors and resolutions obtained in the 
respective frequency components. 

20 Finally, the quantized bases U q § 'Mx), V q * i(y) 

obtained in the quantization processing at step S3 6 are 
entropy- encoded at step S37. In the entropy coding 
processing, all the quantized bases U q f 'S.(x) , V q ' i(y) 
are handled as data to be coded. As the method of 
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entropy coding, well-known coding methods such as the 
Huffman coding or Arithmetic coding may be employed. 
Then, the entropy-coded quantized bases, the number of 
bases by window, and index information necessary for 
5 decoding such as window size, are recorded as a part of 
the coded data, and the process ends. The index 
information may be omitted if a previously- standardized 
index information is used. 

In the third embodiment, in the DCT transform 

10 processing at step S35, discrete cosine transform is 
performed directly on plural discrete one -dimensional 
adaptive base groups , Y^i(y) , however the 

present invention is not limited to this arrangement . 
It may be arranged such that the respective one- 

15 dimensional adaptive base groups X?' 6 x (x) . ^\{y) are 
further divided, then the divided base groups are 
subjected to discrete cosine transform, quantization, 
and entropy coding. 

To decode the coded data obtained by the above - 

20 described method to a decoded image f decode (x,y) , the 
quantized bases U, 5 t i(x) , V^\(y) are generated by 
entropy decoding and inverse-quantized while the index 
information is appropriately used, thereby decoded DCT 
vectors Udecode 1 '^*) . Vdecode^iCy) are reconstructed, and 

25 inverse discrete cosine transform is performed on the 
reconstructed decoded DCT vectors, thereby decoded one- 
dimensional adaptive bases Xaecod© i(x) , Y d ecode ' i(Y) are 
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reconstructed. Next, based on the reconstructed 
decoded one-dlmensional adaptive bases X de coae ' x(*) . 
Ydecode ?,E i(y) . decoded windows W ae oode l ^(x,y) is generated 
as follows. 

5 

{x,y)~ ^X^n (*)- Y^J** (y )] 

In the above expression, k 1 ^ is the number of 
selected and corrected base groups by window. The 
10 decoded image f decode (x,y) is obtained by combining the 
generated decoded windows Wdecode^fx.y) generated as 
above . 

In the third embodiment, a difference coded image 
is generated by obtaining a reconstructed image from 

15 the generated discrete one -dimensional adaptive base 
groups X^i(x), Y ? 'S.(y) as in the case of the second 
embodiment. As the error is absorbed by quantization, 
it may be arranged such that upon generation of each of 
the sequentially generated one -dimensional adaptive 

20 base groups, the DCT transform at step S35 and the 

quantization at step S36 are performed, then inverse 
quantization and inverse discrete cosine transform are 
performed on the quantized bases U q ? * 5 i(x) , V q l '^i(y) by a 
method similar to that of the decoding, thereby the 

25 reconstructed decoded one-dimensional adaptive bases 
Xdecode^i(x) . Ydecode^iCy) are generated, then a 
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reconstructed Image is obtained by using the 
reconstructed decoded one -dimensional adaptive bases, 
thereby a difference coded image is generated. 

Fig. 28 shows an example of a decoded image 
5 obtained by applying the above coding and decoding to 
the image in Fig. 24. In the coding, the image is 
divided into 8x8 pixel blocks (Lx=Ly=8) and coding is 
performed by block. In this example, the DCT transform 
and the quantization are performed upon generation of 
10 each one -dimensional adaptive base group, then inverse 
quantization and inverse discrete cosine transform are 
performed, and a difference coded image is generated by 
using the result of inverse DCT transform. As the 
error due to quantization is absorbed upon generation 
15 of next one -dimensional adaptive base group, the 

accumulation of quantization error can be prevented. 
Regarding remainder (s) which occur upon division into 
8x8 pixel blocks, coding is performed on the assumption 
that the image has all 8x8 pixel blocks and the values 
20 of portions other than the image area are "0". Further, 
the quantization coefficient vectors Q x (x) and Q y (y) 
are as shown in Fig. 29. In this case, the data amount 
is, including information necessary for decoding (image 
size: 4 bytesx2, the number of base groups of each 
25 block; 1 bytexthe number of blocks [13x13 groups], 
quantization coefficient vectors: 4 bytesxS), 3790 
bytes. That is, the data amount is further reduced 
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while the image has quality equal to or higher than 
that of 8 -bit rounded image in the second embodiment. 
In this example, the image is divided into 8x8 pixel 
blocks and coding is performed by block, however, other 
block size such as 20 x 20 pixels may be used upon 
division. 

In the coding method described in the third 
embodiment, as in the case of the second embodiment, an 
input image is divided into plural windows and base 
grouping is performed, thereby the number of base 
groups corresponding to each window can be variably set. 
Accordingly, in a window having monotonous density 
distribution as shown in Fig. 5, coding with reduced 
error can be performed with a reduced number of base 
groups. On the other hand, in a window having density 
distribution of natural image or the like, coding can 
be performed with a large number of base groups. Thus 
the number of base groups can be optimized. Further, 
as the base groups are discrete cosine transformed then 
quantized and entropy encoded, higher efficiency coding 
can be realized. 

As described in the third embodiment, a discrete 
two-dimensional gray-scale image f(x,y) defined in a 
finite range is divided into plural blocks, then plural 
discrete one -dimensional adaptive base groups 
Xdecoae^iU), Y dec ode i ' 5 i(y) are generated by block, and the 
base groups are encoded. 
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< Fourth Embodiment > 

In a fourth embodiment as a modification to the 
third embodiment, a discrete two-dimensional gray-scale 
5 image f (x,y) defined in a finite range is subjected to 
two-dimensional discrete cosine transform and 
quantization, thereby a quantized frequency image 
f u ' v i(?£) is generated, and encoded to plural discrete 
one -dimensional adaptive base groups X u ' v x (x) , Y u ' v i(y) . 

10 That is, the two-dimensional image is converted to 
frequency component data (represented in frequency 
space) and subjected to the above-described coding 
using adaptive bases. 

In the fourth embodiment, the apparatus 

15 construction is the same as that of the first 

embodiment. Hereinbelow, the processing according to 
the fourth embodiment will be described with reference 
to the flowchart of Fig. 4. 

First, in the image input (step S40), an input 

20 image f(x,y) as a discrete two-dimensional gray- scale 
image to be encoded is inputted. In this embodiment, x 
and y have limited definition ranges similar to those 
in the second embodiment, and are integers as in the 
case of the second embodiment. 

25 Next, two-dimensional DCT transform processing is 

performed (step S41). The input image f(x,y) is 
divided into plural images W^(x,y) having a 
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predetermined size of, e.g., 8x8 pixels. The divided 
images W^(x,y) are similar to the windows W^(x,y) 
described in the third embodiment. At step S41, two- 
dimensional discrete cosine transform is performed on 
5 the respective divided images W^(x,y) . 

Next, the two-dimensional discrete -cosine 
transformed data are quantized (step S42). The 
quantization processing is performed by using a 
predetermined quantization coefficient Q(x,y) , thereby 
10 quantized images W q u (x,y) are generated. That is, 
processings similar to the JPEG coding step of two- 
dimensional discrete cosine transform and quantization 
are performed. Since the transform is similar to the 
JPEG transform, the detailed explanation thereof will 

15 be omitted. 

In the fourth embodiment, the input image f (x,y) 
is divided into 8x8 pixel sized images W l,e (x,y) and 
subjected to two-dimensional discrete cosine transform, 
and subjected to quantization by using the quantization 

20 coefficient Q(x,y) to increase a value in a high 
frequency, however, the present invention is not 
particularly limited to this arrangement. For example, 
the input image may be divided by another size, or the 
quantization may be performed by using another 

25 quantization coefficient. 

Next, frequency selection processing is performed 
(step S43). In this processing, frequency parameters u. 
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v are selected, and a frequency image f u ' v (i,£) is 
generated. The frequency parameters u, v are integers 
equal to or greater than "0". As the input image 
f (x,y) is divided into 8x8 pixel sized images W^(x,y) 
5 in the two-dimensional DCT transform (step S41), the 
parameters are integers up to "7". Accordingly, the 
frequency parameters u, v have values "0 M to "7". In 
this embodiment, one of 8x8=64 parameter combinations 
is selected. The generated frequency image f u * v (££) is 
10 expressed as follows using the quantized images 
W q ^(x,y) generated in the quantization. 

f u v (££) - W q ^(x,y) 

15 As expressed above, in the frequency image 

f u ' v (£,C) . selected u, v frequency components in the 
respective quantized images W^x.y) are arrayed. 

Next, in the base grouping processing (step S44), 
the frequency image f u ' v (l£) generated in correspondence 

20 with the selected frequency parameters u, v is 

converted to plural discrete one-dimensional adaptive 
base groups X u '\(i) , Y a -\(t,) . Since the processing is 
similar to the coding to plural one -dimensional 
adaptive base groups described in the second embodiment, 

25 the explanation thereof will be omitted. 

In selection completion determination (step S45), 
it is determined whether or not all the frequency 
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parameters u, v have been selected. In the fourth 
embodiment, as there are 64 combinations of frequency 
parameters u, v, all of the combinations are selected. 
When the base grouping of the frequency images f u * v (|,U 
5 corresponding to the respective frequency parameters u, 
v in the base grouping (step S44) has been completed, 
the process proceeds to entropy coding (step S46) . If 
it is determined that all the 64 combinations of 
frequency parameters u, v have not been selected, the 
10 process returns to the frequency selection (step S43), 
in which an unselected combination of frequency 
parameters u, v is selected, then base grouping of the 
frequency image f u ' v (S£) corresponding to the selected 
frequency parameters u, v is performed. 
15 m the entropy coding (step S46), the base groups 

X u,v i(|). Y u ' v i(£) generated in the base grouping (step 
S44) are entropy encoded, then recorded, with index 
information necessary for reconstruction, as coded data, 
and the process ends. The processing is similar to 
20 that described in the third embodiment. 

As in the case of the third embodiment, in the 
method for generating a decoded image from the coded 
data, first, one -dimensional base groups are 
reconstructed by entropy decoding with the index 
25 information, then a frequency image is reconstructed by 
using the base groups. Then quantized images are 
generated from the frequency images, then inverse- 
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quantization is performed, and inverse two-dimensional 
discrete-cosine transform is performed, thereby a 
decoded image can be reconstructed. 

In the third embodiment, the input image is 
5 spatially divided and encoded, while in the fourth 
embodiment, the input image is divided in frequency 
space and encoded. In the fourth embodiment, one- 
dimensional base grouping is made from the image 
divided in frequency space, however, it may be arranged 

10 such that the image divided in the frequency space is 
further spatially divided before base grouping. 
Further, in the entropy coding (step S46), the 
generated base groups are directly entropy encoded, 
however, it may be arranged such that the one- 

15 dimensional adaptive base groups are quantized and 

entropy encoded. Further, as in the case of the third 
embodiment, it may be arranged such that further 
discrete cosine transform (one-dimensional) then 
quantization are performed, and entropy coding is 

2 0 performed . 

As described above, a discrete two-dimensional 
gray-scale image f (x,y) defined in a finite range is 
two-dimensional discrete-cosine transformed and 
quantized thereby a frequency image f u ' v (§,£) is 

25 generated, and encoded to plural discrete one- 
dimensional adaptive base groups X uv i(£), Y u ' v i(£) . 
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< Fifth Embodiment > 

In the above -described first to fourth 
embodiments, coding of two-dimensional data is 
described. However, the present invention is not 
5 limited to two-dimensional data but is applicable to 
other multidimensional data. In a fifth embodiment as 
a modification to the second embodiment, discrete 
three-dimensional data f(x,y,z) defined in a finite 
range is encoded to plural discrete one -dimensional 
10 adaptive base groups Xi(x) , Y x (y) , Zi(z). 

As discrete three-dimensional data to be encoded, 
discrete density distribution data in three-dimensional 
space or discrete time-series data of discrete two- 
dimensional gray- scale image may be used. In this 
15 embodiment, each element of the three-dimensional data 
(in the case of discrete time- series data of discrete 
two-dimensional gray-scale image, a pixel of the two- 
dimensional gray-scale image at predetermined time) is 
represented by 1 byte. 
20 First, such three-dimensional data f(x,y,z) to be 

encoded is inputted. As in the case of the second 
embodiment, x, y, z are Integers having predetermined 
definition ranges 0sx<x max , 0ssy<Vma*, and 0sz<z max . As in 
the case of the second embodiment, coding is performed 
25 by correcting one -dimensional adaptive base groups 

X x (x) , Yi(y), Zi(z) based on the following error E and 
correction amounts AXi(x) , AY x (y) , AZi(z). 
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In the above expression, fi(x,y,z) and f t (x,y,z) 
5 are similar to the difference coded image fx(x,y) and 
reconstructed image f t (x,y) in the second embodiment. 
The only difference is that the dimension is three- 
dimensional. Since other parameters are similar to 
those in the second embodiment, the explanations 

10 thereof will be omitted. 

The three-dimensional data is encoded to the one- 
dimensional adaptive base groups X x (x) , Y x (y) , Zi(z) by 
correcting the one- dimensional base groups based on the 
error E and the correction amounts AXi(x) , AYi(y) , 

15 AZi(z), by a similar method to that described in the 
second embodiment. Assuming that one -dimensional 
adaptive base group corrected from selected one- 
dimensional adaptive base group Xi(x) , Yz(y) , Z x (z) are 
X'i(x), Y'i(y), Z'i(z), the correction is performed as 

20 follows. 
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z;(z)-z,{z)+*z f (?) 



Finally, when the error E is equal to or a 
5 predetermined value, or the number of bases becomes a 
predetermined number, the one -dimensional adaptive base 
groups Xi(x), Yi(y) , Zi(z) selected and corrected by 
that time are recorded, and the process ends. Since 
the processing here is similar to that described in the 
10 second embodiment, the detailed explanation thereof 
will be omitted. In this case, decoded data 
fdecode(3c,y,z) reconstructed from the one- dimensional 
adaptive base groups X x (x) , Yi(y), Zi(z) is as follows. 

15 CT 

In the above expression, k is the number of base 
groups. In this manner, even if the dimension of data 
to be encoded becomes higher, coding can be realized by 
20 obtaining error and correction amounts as described 
above. In the present embodiment, the input three- 
dimensional data is directly encoded to one -dimensional 
adaptive bases, however, as in the case of the above- 
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described third embodiment, the three-dimensional data 
may be divided by a predetermined range or divided in 
frequency space by discrete cosine transform. 

As described above, discrete three-dimensional 
5 data f(x,y,z) defined in a finite range is encoded to 
plural discrete one -dimensional adaptive base groups 
X!(x), Y x (y), Zi(z). 

< Sixth Embodiment > 

10 In the above -described first to fifth embodiments, 

as an error evaluation function, the second power of 
error of each element is employed, however, the coding 
method of the present invention is not limited to the 
above evaluation function, but is applicable to other 

15 error evaluation functions. In a sixth embodiment as a 
modification to the second embodiment, assuming that a 
is reference data, i.e., data to be encoded, and b is 
evaluation data, i.e., data reconstructed from coded 
data, an error evaluation function 2{ l/2x(a 2 -ab) 2 } is 

20 employed when a discrete two-dimensional gray-scale 
image f (x,y) defined in a finite range is encoded to 
plural discrete one- dimensional adaptive base groups 
Xx(x), Y^y). 

First, a discrete two-dimensional gray-scale 

25 image f (x,y) to be encoded is inputted. As in the case 
of the second embodiment, x, y are integers having 
predetermined definition ranges 0:sx<x ma x# 0^y<y m ax- With 
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regard to the two-dimensional gray-scale image f(x,y), 
coding is performed by correcting the one -dimensional 
adaptive base groups Xi(x) , Yj.(y) based on the 
following error E and correction amounts AXi(x) , AYi(y) , 
5 as in the case of the second embodiment . 



10 similar to the difference coded image fi(x,y) and 

reconstructed image ft(x,y) in the second embodiment, 
therefore the explanations thereof will be omitted. 
Further, since other parameters are similar to those in 
the second embodiment, the explanations thereof will be 

15 omitted. 

Based on the above error E and correction amounts 
AXi(x), AYi(y), the one -dimensional adaptive base 
groups are corrected in a similar method to that 
described in the second embodiment , thereby coding to 
20 the one-dimensional adaptive base groups X x (x), Y x (y) 
is performed. Assuming that one -dimensional adaptive 
base group, corrected from selected one -dimensional 




In tne above expression, fi(x,y) and f,(x,y) are 
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adaptive base group X x (x) , Yi(y) , is X' i(x) , Y'i(y), and 
uncorrected base group is Xi(x), Y x (y) , correction is 
performed as follows. 

AAf,(*) 
5 ^ / (y)»^(y)+Ay / (y) 

Finally, when the error E is equal to or a 
predetermined value, or the number of bases becomes a 
predetermined number, the one -dimensional adaptive base 
groups Xi(x), Y x (y) . selected and corrected by that 
time are recorded, and the process ends. Since the 
processing here is similar to that described in the 
second embodiment, the detailed explanation thereof 
will be omitted. In this case, decoded data f decode (x,y) 
is similar to that described in the second embodiment. 

In this manner, coding is realized by using other 
error evaluation function than the second power error, 
by obtaining error and correction amounts as described 
above. In the sixth embodiment, the error evaluation 
function 2{l/2x(a 2 -ab) 2 } is employed, however, the 
present invention is not limited to this error 
evaluation function. Any other evaluation function may 
be used as long as it becomes closer to coded data 
thereby error is reduced. For example, in the case of 
two-dimensional image, the divergence between 
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difference coded image and differentiated reconstructed 
image may be included in the error evaluation function. 
It is preferable that in such error evaluation function, 
parameters to be corrected, i.e., partial differential 
5 with respect to Xl(x) , Yi(y) is analytically obtained. 
However, if the partial differential cannot be 
analytically obtained, the correction amounts may be 
determined from increment /decrement of error when the 
parameters to be corrected are slightly changed. 

10 As described above, a discrete two-dimensional 

gray- scale image f(x,y) defined in a finite range is 
encoded into plural discrete one -dimensional adaptive 
base groups Xi(x) , Y x (y) by using 2{l/2x(a 2 -ab) 2 } as an 
error evaluation function. Note that the sixth 

15 embodiment is particularly applied to the second 

embodiment, however, it is also applicable to other 
embodiments . 

< Seventh Embodiment > 

20 In the above-described first to sixth embodiments, 

multidimensional data is encoded to plural one- 
dimensional adaptive base groups, however, the present 
invention is not limited to the above coding, but is 
applicable to coding to other form of data. In a 

25 seventh embodiment as a modification to the second 

embodiment, a discrete two-dimensional gray- scale image 
f(x,y) defined in a finite range is encoded to plural 
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discrete two-dimensional adaptive bases S x (x s ,y s ). 

First, a discrete two-dimensional gray-scale 
image f (x.y) to be encoded is inputted. As in the case 
of the second embodiment, x and y are integers having 
5 predetermined definition ranges 0 sac<x max and 0<y<yma*. 
Next, as a base, one discrete two-dimensional adaptive 
base Si(x s ,y s ) is selected. If x,^ or ymax is an even 
number, x s and y s are integers having definition ranges 
0sx s <(x ma x/2)+l and Osy s < (yma*/2 ) +1 . If Xmax or y max is an 

10 odd number, x s and y s are integers having definition 

ranges 0*x s < (Xmax/2) +0 . 5 and 0^y s < (y»ax/2 )+0 . 5 . Further, 
the suffix "1" is a base number. Although not 
particularly limited, the base number is sequentially 
selected from "0" in the present embodiment. 

15 Next, a reconstructed image f,(x,y) is calculated 

from the selected two-dimensional adaptive base 
Si(x s ,y s ) having the base number "1" as follows. 
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i ) y : even number 

ii ) x s y : odd number 

lii) x : even number, y : odd number 

iv) x : odd number, ? = even number 

The error E and correction amount ASj.(x g ,y s ) are 
calculated by using the above reconstructed image 
5 f t (x,y) as follows. 

In the above expression, fx(x,y) is a difference 
10 coded image similar to that of the second embodiment. 
Further, in the expression to calculate the correction 
amount ASx(x s ,y s ) , C(m,n) has the following values. 
If n=0 and m=0 hold 

C(n,m)=l 
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If n=±l and m=0 hold, or n=0 and m±l hold 

C(n,m)=0.5 
If n=±l and m=±10 hold 

C(n,m)=0.25 

5 

Note that in the above expressions, if 2x s +n<0, 
2x s +na3Cmax, 2y s +m<0, or 2y s +n2 : y ma2: holds, C(m,n)=0 holds. 

Further, r\ is a correction parameter similar to 
the parameters y\ x and ■% described in the second 
10 embodiment. The correction is performed on the two- 
dimensional adaptive base based on the above error E 
and correction amount ASi(x s ,y a ) in a similar method to 
that described in the second embodiment, thereby the 
input image f (x,y) is encoded to the two-dimensional 

15 adaptive base Si(x s ,y s )- 

Regarding the correction of the base, the 
selected two-dimensional adaptive base Si(x s ,y s ) is 
corrected as follows. In the expression, Si(x s ,y s ) is 
an uncorrected base prior to the correction, and 

20 S'i(x s ,y s ) is a corrected base. 
S , i (x s ,y 3 )=Si(Xs,ys)+ASi(x a ,y s ) 

Finally, when the error E is equal to or a 
predetermined value, or the number of bases becomes a 
25 predetermined number, the two-dimensional adaptive 

bases Si(x s ,y s ) selected and corrected by that time are 
recorded, and the process ends. Since the processing 
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here is similar to that described in the second 
embodiment, the detailed explanation thereof will be 
omitted. In this case, a decoded image f decode (x,y) 
reconstructed from the coded data is obtained as 
5 follows . 

i ) y l even number 

ii ) y : odd number 

ill) x : even number , y : odd number 

/^>)-l[i{^)^i.^ 

i\r) x : odd number, y : even number 

In the above expressions, k is a base number of 
10 selected and corrected bases. In the seventh 

embodiment, in the generated two-dimensional adaptive 
bases S x (x s ,y s ), the number of elements is about 1/4 of 
the input image f(x,y), however, the invention is not 
limited to such number but it may be arranged such that 
15 coding is performed by using a two-dimensional adaptive 
base having a different number of elements. 
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As described above, in the seventh embodiment, a 
discrete two-dimensional gray- scale image f(x,y) 
defined in a finite range is encoded to plural discrete 
two-dimensional adaptive bases Si(x s ,y s ). 

5 * 
< Eighth Embodiment > 

In the above-described first to seventh 
embodiments, coding is performed by using a one- 
dimensional adaptive base or two-dimensional adaptive 
10 base as an adaptive base, however, the coding method of 
the present invention is not limited to this 
arrangement but is applicable to adaptive base having 
plural parameters. In an eighth embodiment of the 
present invention as a modification to the second 
15 embodiment, a discrete two-dimensional gray- scale image 
f (x,y) defined in a finite range is encoded to an 
adaptive base having plural parameters. 

Hereinbelow, as a particular example, coding to 
an adaptive base Gi(Ai,pi,qi,pi) having four parameters 
20 Ai, pi, qi and pi will be described. 

First, a discrete two-dimensional gray-scale 
image f (x,y) to be encoded is inputted. As in the case 
of the second embodiment, x and y are integers having 
predetermined definition ranges 0sx<x max , 0sy<ymax- Next, 
25 one adaptive base Gi is selected as an initial adaptive 
base. The suffix "1" is a base number. In the eighth 
embodiment, although not particularly limited, the 
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adaptive base is sequentially selected from "0" to a 
predetermined number. 

Next, a reconstructed image f~i(x,y) is 
calculated from the selected adaptive base d having 
5 the base number 1 as follows. 

Next, the error E and correction amounts AA i# Api, 
10 Aqi, and Ap x are calculated as follows by using the 
above reconstructed image f"i(x,y) . 

*-?$[£faw>-«wfr] 



6A, = 


BE 

~BA,~ 




APi = 


BE 
-n ■ 

' 9P, 






BE 




Aa, = 


BE 
" Bo, 





15 m the above expression, fi(x,y) is a difference 

coded image similar to that of the second embodiment. 
Further, T] A . Tip , Tj g and t] 0 are correction parameters 
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similar to the parameters Y| x and t] y described in the 
second embodiment. Based on the above error E and the 
correction amounts AA X , Api, Aq if and Api, the 
reS p ec tive parameters constructing the adaptive base 
are corrected by a similar method to that of the second 
embodiment , thereby coding to the adaptive base d is 
performed. 

Regarding the correction of the adaptive base, 
the respective parameters Ai, pi, qi and pi, 
constructing the selected base G x are corrected as 
follows. A f i, p'i, q'i and p f i are corrected parameters, 
and Ai, pi, qi and pi are uncorrected parameters. 
A'i = Ai + AAi 
P f i - Pi + Api 
q'i = qi + Aqi 
p f i = pi + Api 

Finally, when the error E is equal to or less 
than a predetermined value, or the base number becomes 
a predetermined number, the respective parameters A x , 
pi, q x and pi constructing the adaptive bases G x 
selected and corrected by that time are record- 
output ted, and the process ends. Since the processing 
here is similar to that described in the second 
embodiment, the detailed explanation thereof will be 
omitted. In this case, a decoded image f decode (x,y) 
reconstructed from the coded data is obtained as 
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10 



15 



follows . 



f decoded) 



k-1 



<.-f~S£=» t }] 



In the above expression, k is a base number of 
selected and corrected bases. In the eighth embodiment, 
the adaptive base Gi having four parameters is used, 
however, the present invention is not limited to this 
number of parameters , but it may be arranged such that 
coding is performed by using an adaptive base having a 
different number of parameters. 

As described above, a discrete two-dimensional 
gray- scale image f (x,y) defined in a finite range is 
encoded to adaptive base d having plural parameters. 



<Ninth Embodiment > 

In a ninth embodiment, discrete convolution 
operation to two-dimensional distribution data as 
predetermined two-dimensional weight distribution data 

20 is performed. The two-dimensional weight distribution 
data is transformed by the data coding according to the 
second embodiment to one -dimensional adaptive base 
group, and discrete convolution operation is performed 
on the two-dimensional distribution data by using the 

25 one -dimensional adaptive base group. The discrete 
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convolution can be used in, e.g., calculation in one 
feature detection cell face in Convolutional Neural 
Networks ( LeCun , Y . and Bengio , Y . , 1995, 
"Convolutional Networks for Images Speech, and Time 
5 Series" in Handbook of Brain Theory and Neural Networks 
(M. Arbib, Ed.), MIT Press, pp. 255-258). In this case, 
the two-dimensional weight distribution data can be 
regarded as synapse weight distribution. 

First, discrete two-dimensional weight 

10 distribution data w(x,y) to be subjected to discrete 

convolution is inputted. As the two-dimensional weight 
distribution data is discrete, x and y are integers. 
Note that if the two-dimensional weight distribution 
data w(x,y) is k-dimensional weight distribution data, 

15 the data is expressed as w(x 1# x 2 , . x k ) . Further, 

if the two-dimensional distribution data is k- 
dimensional distribution data, the data is expressed as 
A(xi, x 2/ x k ). In the present embodiment, the two- 

dimensional distribution data is expressed as f (x,y) . 

20 The result of discrete convolution to the two- 
dimensional distribution data f (x,y) as two-dimensional 
weight distribution data is expressed as follows. 

25 

Next, a one-dimensional adaptive base group Xi(x) , 
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Yx(y) is generated from the Inputted two-dimensional 
weight distribution data w(x,y) by a similar method to 
that described in the second embodiment. The suffix 
"1" of the one-dimensional adaptive base group Xi(x) , 
5 Yi(y) is a base number. As detailed explanations here 
are similar to those in the second embodiment , the 
explanations will be omitted. The one -dimensional 
adaptive base group is not necessarily generated upon 
each convolution , but it may be arranged such that the 
10 one -dimensional adaptive base group is previously 
generated by the method described in the second 
embodiment, and the generated one -dimensional adaptive 
base group is recorded. 



15 dimensional distribution data f (x,y) by using the 

generated one -dimensional adaptive base group is as 
follows . 



In conventional art as shown in Pig. 6, in 
discrete convolution of two-dimensional weight 
distribution data having 6x6 elements, two-dimensional 
weight distribution data values are independently given 
25 to the total thirty- six elements in the respective 
coordinates, i.e., the two-dimensional weight 



The result of discrete convolution to the two- 




20 
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distribution data values are sent as signals to the 
respective coordinates, thirty- six signal lines are 
required. On the other hand, in the present embodiment, 
as the two-dimensional weight distribution data w(x,y) 
5 is decomposed to the one -dimensional adaptive base 
group Xi(x), Y x (y) and discrete convolution is 
performed by using the one -dimensional adaptive base 
group, the required number of signal lines is twelve as 
shown in Fig. 7. In this manner, the conventional 
10 problem of wiring in the convolution arithmetic unit 
can be solved. Note that in a convolution arithmetic 
unit having an array of product -sum arithmetic elements 
for execution of the above discrete convolution, the 
computation may be performed by inputting the one- 
15 dimensional adaptive base group Xi(x), Y x (y) into all 
the columns and rows of the array, as shown in Fig. 7, 
otherwise, by inputting the one -dimensional adaptive 
base group X x (x) , Y^y) into respectively only a part 
of the columns and rows, and by performing computation 
20 in only a part of area (63) of the array where both of 
Xi(x) , Yi(y) are inputted, as shown in Fig. 9. 

As described above, in discrete convolution to 
two-dimensional distribution data as predetermined two- 
dimensional weight distribution data, the two- 
25 dimensional weight distribution data is converted to 

one- dimensional adaptive base group by the data coding 
method according to the second embodiment, and discrete 
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convolution is performed on the two-dimensional 
distribution data by using the one -dimensional adaptive 
base group. 

In the above embodiment, k- dimensional data w(xi, 
5 x 2 , .... x k ) is encoded by the above -described coding 
method, thereby at least a set of k one -dimensional 
adaptive base groups Xn(xi) , x 2 i(x 2 ) , x k i(x k ) is 

obtained. Further, the result y of the product -sum 
operation using the above k-dimensional data A(x x , x 2 , 
10 . . . , x k ) and the above one -dimensional adaptive base 
group is calculated as follows; 

15 By the above calculation, the result y can be 

efficiently obtained by product -sum operation between 
the k-dimensional data w(x 1# x 2 , x k ) and the k- 

dimensional data A(x x , x 2 x k ) as follows. 

2 o y =JJ-J^,a: 2 ,...,^)-w(^,a: 2 ,...,jcJ^A 2 ..^ 

In addition to the embodiments of the present 
invention as described above, as it is apparent that 
the present invention can be realized by a computer 
25 program as shown in Fig. 8, such computer program is 
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Included ±n the scope of the present Invention. 
Generally, a computer program is stored in a computer- 
readable storage medium such as a CD-ROM, and when the 
medium is set in a storage medium drive device of a 
5 computer then the program is copied or installed to the 
system, the program becomes executable. Accordingly, 
such computer-readable storage medium is also included 
in the scope of the present invention. 

As described above, according to the embodiments, 

10 n data to be encoded can be encoded to data of 

efficient representation format (format using adaptive 
base). For this purpose, the adaptive base is set, 
then n data are reconstructed by using the adaptive 
base, and the error between the n data to be encoded 

15 and the reconstructed n data is calculated by using a 
predetermined error evaluation function, and the 
adaptive base is corrected based on the error. 

As described above, according to the present 
invention, high-ef f icient data coding, more 

20 particularly, coding of distribution data in 

multidimensional space such as a two-dimensional image 
can be realized. Further, as two-dimensional data can 
be decomposed to one-dimensional base data, the 
conventional problem in convolution arithmetic unit, 

25 i.e., the problem of wiring of two-dimensionally 
parallelized array, can be solved. 

As many apparently widely different embodiments 
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of the present invention can be made without departing 
from the spirit and scope thereof, it is to be 
understood that the invention is not limited to the 
specific embodiments thereof except as defined in the 
5 appended claims. 

CLAIM OF PRIORITY 

This application claims priorities from Japanese 
Patent Application No. 2004-024485 filed on January 30, 
10 2004 and Japanese Patent Application No. 2004-239792 

filed on August 19, 2004, the entire contents, of which 
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